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(54) HIGH SPEED TECHNICAL CALCULATION METHOD 

(57)Abstract: 

PURPOSE: To perform a technical calculation at a high speed without 
using any expensive numerical arithmetic processor and at the same time 
to reduce the necessary capacity of a memory byperforming a numerical 
operation through an integer arithmetic and in a block floating point 
format. 

CONSTITUTION: A fixed floating point arithmetic using an integer 
arithmetic is expressed in a block floating point number that showed the 
data in an exponent part or a coefficient part common to all data and a 
fixed decimal point of each data. Then a block floating point arithmetic 
applying an integer arithmetic is carried out. It is supposed that a decimal, 
point is set at the 13th bit of an integer and then that 8192=213 of an 
integer is equivalent to 1.01 of a fixed decimal point. In this fixed decimal 
point format, the value covering -4.0 (minimum) through about 3.9999 
(maximum) can be expressed. The expressible number most approximate 
to 0 is about 0.000122 and a number approximate to 1.0 has the decimal 
accuracy equivalent to about 4 digits. The diagram shows a program used 
for calculation of a complex number end applying a fixed decimal point in 
a C language. 



LEGAL STATUS 

[Date of request for examination] 

[Date of sending the examiner's decision of rejection] 

[Kind of final disposal of application other than the 
examiners decision of rejection or application 
converted registration] 

[Date of final disposal for application] 

[Patent number] 

[Date of registration] 

[Number of appeal against examiners decision of 
rejection] 

[Date of requesting appeal against examiners decision 
of rejection] 

[Date of extinction of right] 

Copyright (C); 1 998,2003 Japan Patent Office 




I 



http://www1 9.ipdl.jpo.go jp/PA1 /result/detail/main/wAAA3IaiqADA4041 77461 P1 .htm 



2004/06/18 



®&m¥fWf2:m (a) ¥4-177461 



©Int. CI. 5 Jrf*)®i#^ ®£§B ^4^(1992) 6 ^ 24 B 

G 06 F 15/31 Z 6798- 5L 

15/332 A 6798- 5 L 

gaging ^n* m*m<»& i (^iqm) 

©%f m ¥2- 305035 

@* li ¥2(1990)11^7 0 



I i 
2 . #atw#<Dtt® 

- * **&JI+*ft*ffit*Ci3V*T. «E « » » €: f«J 

«$ 5 IS ( a ) , (b)tt, CQa5JK<0 r CJC«fc 
D y ^ 

(1) iiajRKfca#**J&5£cosplex££|| 



doubled iaT L T V> * . (2) I* (1) T* 

fell I, feco BP lex^O^»«0B8»caddT«> V , 
(3) i*complexg?<9il*i*flI 0>l»ttcsuL (4) 
coup \e%m CO # <7)|»«[ caul T* y , ( 5 ) 
t* complex^ ft nm O 08 KcdivT? & * . (6) 
ticooplex^ <0 ft 6D*a*f & # 4b £> 08 Efc cabs T & 
(7) tidouble3$0>tt 2 o r * 
coapiex^ O&$::£#>£>0Bft to-coaplexT' tb £> . 

#6|giilS3D< rciC * **43*8cffilt#J 0)21 
8 * - V \ C ft IS $ n T W% ttt * fit » * & 8 V > tz F F T 
^D^5iA-e«>45. (1 0) ttt^^-f 

(i2)«M^7^>fatffiio5fesa«r^ 

#570 (a), (b) tt»60^7ti:tT06 
y^/^T^IS^tlTV^fc FFT« a&<Dm>) 

(i o) - (i 2) ttKeakwct^T 



-395- 



$> *J . ( 1 3 ) H FFTfc »^0ffi^JZ©SJ8$: n 
T fctf & »l k J # 80 £ ^ L/ T V* 6 . 31 8 0 tt . tT y h 
ij^-f/l/Hfl)7D-f ir-K IS 9 ® II * * 7 
7-{«il(07D-f y-K 85 i 0 13 li £i »J * $ 

£) "7 D — f- Y — 

3&5g|0> (2) OKfi[cadd«^58«S[a)See: a , 

Z.re(Z©*»)«». re(a<0«»)*b . re(bCO^gP) 
Z. io(Z<D£a6iO° a . iB(atf)&g{J)+ b . i«(btf>&3fl) 

£ jg f ^ <0 T . X. re* a . re- b . re, X. im = a . io- 
b.iofeft*LT^&. 

£15 12 03 (4) OflBKcounil^^ttCUl^SaCDSl 

X. re= a . reX b.re- a.ioX b.io 
X. io= a . reX b . in- a.ioX b.re 
* if-* U v**. 



J5fin 2 M-177>lGl(2) 
§5500 (5) ttWBcdivttra8UC#$&<oi$ 

a. reX b . re + a.ioX b.im 
X. re= — 

b. reXb.re+b.inXb.lo 

a. imX b . re + a.reX b.io 
X. in- ~ 

b. rex b.re+ b.ioX b.im 

£ It* LTV* 6. X, #S#0<£>B*tix5>-*i* 

f&SEItf) (6) CBBB&cabstt a <D*S*J"fI £ double 
T'ig1-M#-e I a | W( a . re) * + ( a . io) * 
If-* LT . 

fl 5 I2tf> (7) 0>MBt tocooplex«doubleC03CK 
ICSCtft t Tig* B8»T Z . re«= x , Z.in»y*% 
«6gl, * 7 Htttt3R»W*fe«^fcFFT©^ 
#8l2J*itJfia<Z>:/ny^i*tf>try hU^-t^ 



meiattiiasB^p^^-kcD^* ? ^ -f ?» * o> 
7D-ft-ht»y, m w<d (i 2) <*> an # 

$100(^700 (13) (#6l3T'U:*tt 
S ft T V> * ) (0|!»Jfli|ro7D-f^-hTJ> 

K^JZ. FFTCO^IR 1 . FFT il BFT £ 

aHR***»<&3CKf ( f - 1 * 6RFT, f«- 
1 * £FFT#Htfa ft 6) HIRtU^iT^ 
*7 f « 1 e»fi2^i Z Ot 5 - * 7 - "J 

- I fx t> s fi^JZOf-^ fc7-'JxXlftlfcf 
-*#82?UZU::£*b*>ft*. 

FFT^D^^^H. JUS* 7 ~ U xSCm £ ifcfcT 

g (k) <D^®»^^^h;i/G (n) £ , 

G ( n ) « - £ g ( k ) Wh 
N k -° 

£v*e>J»<&£/8v*T#a*yD^ATifca. <H 



Ntt^-JBT&SZCO^^&Tfcfcaft* 
ft (2,4, 8, •••,128.256, 512, 1024, •••) £. L ♦ 

frftTv*fc*ofc#. ^vn^i^iftJfe 
mSEKD (2) - (7) 0>fc5RB«* £&J8 L, 

rarjSKFFT&K^-r^ i:, m 7 me>± e> kjsv* 

S7Ea(Dya^7At/tV3> PC9801RX2±T 
tK- ? y KttOC 3 > A-f ^ Turbo-C£ flj V* T ^ 
*TSit6fc, X - * » N » 1024CO^r. fijiflia*^ 

»^«* y D -fe v * 180287-10$: 6 t , ISff 
U»*^JW3^THfTTS6, b 180287-1 Oli 
^«T(i 4 77R^^6^P3^^t«<ffi-e*>6. 



-396- 



£64fcr^h<Dttft:&*-i?&9iL , tv*e;&, 10 

iiftTii 1 5rifcOtftt&##6*v* # * 

UK »«5^*a>»*fa*iK»T? *»ft«*<Z> 

IOO~1OOO{&0> #03*51 Ulf-#T6^Wh0 

#>IZ fc 3 ttfc fctf>T, #fffifcftteifc?*:7D iz v * 
O & £ * y©«a%i^<ieV^WieaWlt»y 

Qg/hft A * ^ ti zr a y 9 » a 'J* ft A 



?$nil¥4-17746l(3) 

fiftiHI* £*JJ8 UfcK£/J>ftAi»f* L < tt. ^ 
□ v 0 m tt /J* « & T 8f * £ « & -5 tf) T & * • 
(#*!) 

\z * v^ M £ £ -t * fflr ttf S # 'h ft A & £ <o t> 
* r £ <Z>®5£/i>ftA*e3fe*> L -/D -y * #fiM^ft 

/hKAWW-eff * jikTansimtKaftu* 1 

[JtffcfU 

& & # • & X. \t /t - y 3 > ^ a - * PC980 1 £/ 
'J-XffidDCa>/tO Turbo-C ver2. 0 ( tK - 9 

y K tt ) t ti » fift tt i 6 yd v h a> & # «$ ft 



y hOfff^ttgfiftTl^^jSyofttf^Si-C* 
SMn-32768~* *32767#&3l V £ 5 . 

* ^ & ey r » r tt ft (o i 3 tr v h s k /j> ft a 

tf&*fc<D£<K£UT»ft<0 8l92« 2 ia #ffi£/J> 

ft * coi.oiicra^t-* *aa& ^ fc. ^coa^ 

'J> ft A 7 * - T y h T ti * « /h -4. 0~ « ;* W 3. 99 
99tf 3* #3i-e&*0lC*fci£^ftli#J 
0. 00 01 22 T U 1. O#j5 0> ft <c €j 1 O xg T'W 4 «i 

$> ^ » m tfi th * . 

x^S^/jNftATOiKUtt, ^«m^6fifft2: 

±<wi:*infjieT**^* mti*i.oxi.o 

•Ctt, ^»ftfe32^^h«P^#fi[»^ 
long int*:!*-*-) tCifiC 8t»-T *^5I/)^^. 
3 2 tr^ h^^^«»Tti»/h**J-2«-«^^J 
2 1 ft* -C<Bft * *. long intT^»^ 

8 1 92 T m »J . intl»:g tt ^nif^^^t ^ 



^. «Atf. S«/hft Aft a = 2.0, b =0.5O« 
$:@^/J^ftA<7>CtC^Ar^lS»«:CmgS'e^ia 

Int B,b,c; 

a=(int)(2. 0*8192.0); 
b«(int) (0.5*8192. 0); 

c«(int)((long lnt)a*(long lnt)b/8192); 
fc^e. m b (int)tt int^<7>S[3Clft $: ^ ^ ^ & * 
Xdong int)»ilong int^ 

S: ffi^-T ^ -V ^ h « » ^ T' i . 
/Jn»A«)«I 'J KfeHflllClong intTltflt-T 

ffittfa -1.0, b-2,0Tag4bic 

ic ft a -r 6 li 

int a ,b,c; 

a- (lot) (1. 0*8192. 0) ;b« (int) (2. 0*8192. 0) ; 
c»(int) ((long int)a*8192/(long int)b); 

-(^i:^^®^/hftA^ffl»:€:«ft»»^«a^ 
$iblC»5R»£2oO@^/jN ft Aft tr ^^1*6 



-397- 



r 



rfcfc*rfliT&*. o*yajR#coSS8P&io<o 
Sfc/hB&BtT, 5 I o cog 

t^^tezu^n. c-ffiSTbcoop 

typedef 8truct(int x; int y;)bcoop; 
bcoap Zl ,Z2,Z3; 

x(o»^tt3Rttz i as # z i . x ic . & ap 

Z 1 . y K*ftf-f *. Z1-0.7+J0.7, Z2 
= 0.5+j0.8CO*afcZ3*CttA-f&fC«:±CO£fc 

Zl.x«(lnt)(0.7«8l92.0): 
Zl. y«(int) (0. 7*8182. 0) ; 
Z2.z-(iot)<0. 5*8192.0); 
Z2. y«(int) (0. 8*8192. 0) ; 
Z3. x«Zl. x+Z2. x; 
Z3. y-Zl. y+Z2. y; 

. s^»cogt&£*co**&i£'r*i:*xj!,tc 



?$fin¥4-1774Gl(4) 

*$ i ® (a> , (b) licmaa-csfc/httA** 
co&m8Hi±fc&^cofccofc£<HC?&s. 

(2) COS8»bcaddtt^*/acofl8»-e;fc »J . ft* 
CD & ffi $ <0 B8 « caddtO coaplexCO SP # «: bcoapC ^ 

a <o t & * • 

(3) <OM»bcsubtti*#/8<Ol»»T& U . X 
ftfct*3fc<0&&^c058Rc3ub(OcoaplexCOffil#£ 
bcomptC SlAfc^OT**. 

(4) CO 58 R beau Hi H # * B8RT th V • 
t£ 3fe CO II SSfll CO 58 RcaulCOcoapl exCO fi&# £ bcoap 
ICSCA, tfffl>&8PC0f&* 
22TfifR<0£! €: int/P 6 long i n UC 3£ & b T # •=> 

V>®£>hR,&tC£*&* 8I92Ti0oT 

(5) coBBRbcdivliifc £ # tf>*5SR-e& V * 
t£ & CO 01 CO 58 Red! vCO coaplex £ bcoap IZ & 

?2Rco£3e&&f*fcoTv**. ®£/hR&co 
&T&S<DT«J»»'\8192&&fcT {$ R 7? 



KCO«t<03m-C*y, fli*.** (2.0+jO.O)/(2.0+J 0. 

o > & it * -r * » £ # * <o m m \t n r -c * w t* ( 1 e 3 

84) a + 0«268435456T?& U . *l IC 8192 & & T 5 
#12. 2X 1 0 1 * U long IntTMi^^T £ < 
6. long int<0«*fcttiW2X 10'Tifc »J . % CO fll 
CO 268435456CO 8 f& T S> <& CO T? ^ CO 9t J» T t 
BBfcR2fctf BfcR<0*fi#ia# 2 «Tfc V*}& 

Kttzv* #*ti8te*-etfBrfcT&*. f: 

T, 8«bT^»«8192-*-8-l024-e«|oT^ 
»b«iy#fc&frL. f ©«long in t*> «b In i3? ^ 
lti:^i:«XOfc BfcRco*fittiS#8 

«£TV r n572f«32S^/jNtt^T0. 004i: U /J> gV^I 

^ - £ * * . z <o t$ , x7-*^&t6ie>i: 

#T#fc£8i92-e3UtW*&v%#, XCO^ttiSB 
COJtetttetfftRT / 8192*5 9 0 S £ /J > R A T 0 . 0 1 <fc 
»J /J> $ v> 4: x 9 - t fe 6 . 



2 » l + (8182) a = 32<Oifc*SK LT^^& 3 2 f& L , 
«:256T?« o T * €>S0 y H Sf7«:^.«av^ 
- CD® ^ «?»COJft^<t^H»:T/"256«16, 

/MR^T?0,002J; y /J\$ v^x^ - i: «c ^ . 

l&lEMCfcwT (6) IIIIRttcOie^ti 
tt> 6 88»bcaba-e*<5. t« 3fc <0 £ 16 fll <0 cabscO con 
-plexCO fiP ^ ^ bcoapC fe^3C^/c^COT**>*. 

ff20ICfi^t (7) tt, te^«/^2J^TCO 
dooblecoBC 2 O «: ^n^tl^fiO ,A85 t « bcoap 
(0»fe#^658«ttobcoapTfey , t«*cO|l^05 
CO tocoaplex^: C7?fe6^, ■tt/MtA7 

h Clft'r^C8192«:»ttT v**. 

(8) ttfc JR»bcoapcO»COK^J Z ( ) <0 n 
m<Z>mm& 2TS0e^ >bcnoraT«> U . 

(9) ii^ C < bcoapcoE#J Z ( 3 CO n ® l 



-398- 



aiCZ ( ) <D ±m& £ £8 cons tT «J *J . E 
2 C ) i3D^ji<0|8ttS5lC4S^r 63CBa«p€:con 

J:»J<fc^D^ » »/M»* 7 * - v * h I: it # 
(10) B:, ±TiE^t«URW»«:«mUActt 

*r ft o x v> tfi % a # ^ i*T-tt/i * y ^ << m 

bcnon Z&iT2-em&Z£lZ**)+ N T <D S3 »J 

^ iffitffte>ftic^-*(&**it#W2f8 
/<-:7D-r*rfc#ftv>*^lcfc[i]2T?#j6;i 



«Iin¥4-177461 (5) 

®1 2 0 IC ^ -T . »i 1 S3 JC *3 v* T (1 3) 
tlxlifll $ ft bcnorm-^ yjU-^^^IMfCfcilfca* 

^ J±^M o r v* a . 

69 1 2 I* bcnore-tr ^/b-f >C57*^U X A £ 

9 T v\ 6 . z <n U iz m & /h ft A 7 * - v *y h <n # 

»4Btt. ±I2<&aj»RFFTy O y^Afc-r* 

h-tsaitc^^ftfecmaio^-fvyny^^T 
fey, 10240 0><r - * n&to&y- 4 > * - y £ *s 

-*£A;frLfc8fr<&«ifcR;*'<* h7^i#«>^ 



/=.IS& 1024Atf>FFT#tfj 2«>T^7lfe. 

fr**JRfcft*. X»»t»i(400 K» h X640 

-i/3>^D^HT*;, r ft £ <*> » ft * ft 
+ # T * * . 

SB ft L fc # PASCALS <D flfi tf> ** £ ft L tz M tft ^ SI T 

X, ±i2&»WTtiFFTfcH*7*a:/n^A 
& K ft L ft tfi , mc * $ *Ql&&<»'r- * Z&m 

i* 6 a ffi it » ^ a ^^Afte>Pfl4»<D^ffiTwaHfc 
/>< t£ T & & • 



fttt*J8v%TIt*i-*xfc*;^«BT;fc6. 

X. ±8B^tfc0gT&mUfc-r-*0>**#R* 
3i, tt*5Ri«#&E v*T^«0&ffiIt*£ ! l7ft ^ 
X £ oJtBTJfe 6 a W flt » CO M » ^ o * « 

«*4iC*eft»iSfe*C*»'&^*6«>Ttt* 

&±tf> J: d K £ <Q$fif¥nc Jifttf. Rtem*£fi 
R»Wfc«JBl/T:/D *y *#tt/J>RA? * - ^*y 
hTMiftofcOT, WflSftRtflW»^Dir^y-9-S: 
ft^:tii<, *aT?a«tt*fcfift 3 - * 

4 . ®®<DfW^ft 

®iE!»ira)»^(0-^ifi«JCi«)ia^R^^ 

«*»t»^«fe#«)*iaR«>*ttfea%"ria, * 



-399- 



m 7 muvt&a&m ttfct*8»FFTig. # s em* 
ts 3ft & z <z> » m <?> - m m m & a <» a n » f f t y 

D y 5 2* 0> fcT «y h'J^-*^B©7D-f 

»90«^3R^tt^RFFT^ay^A<Z)^^ 

h. m i i w\t z. o> & m v> - * ^ tft « » 

FFT :/ n^^^CDi^^^-^V BCNORHtf) 7D-f 
Y - h T 6 , 

(i) i*t*3B»0>* (2) 
<fl an J* E IS « . (3) tt1«i&tt<Z)Bt*EW»u 
(4) l*3*3Rtt<Df&*«BBSU (5) 

ass # m m » . (6) u&&&<o&ttmtc##>&w 



?SR1¥4-177461(6) 

(8) tt8 3RBC0>H2?iJ(& n iOS* 2 T? SJ 6 * 
(9) »it$[58jBrOffi^J^S*<0^^ 

*> 0) « * ffi £ 1 <fc »J >J> 3 < ** & * ^ * - * > • 
(10) JR 8 FFT-* y ;i/ - ^ y . ( I 1 ) I* 
fcf*hUJt-*Jl*»* (12) B/t*7M> 

sik (i3) tt«y»»fe»t. 

ra — «f ^ ra — * u < ti*§3&# 

ftiA * £ # m 



T)-in*daf 


nrvct ( 






tat t; 






im r; 












bc**Jt.. b> 




Kuf 


a. b; 












c. « • a. I ♦ b. m 






cr-a.r»b,i; 






r«i«r»( c ); 






faoubta, b) 




boo* 


«. b: 






booap c; 






ci*ci-b.i: 






c. f ■ ». r - b. t: 






miure ( e ); 






bcwtU. b) 






a. b: 






t*xmc c; 






loaf u. t?; 


(ioBc)b.r: 




u • <loai)a.i • (iau)b,i - (loot)*.; • 




t t - (loot)*, i ■ OaM)Kir • Uce«)«.F • 


(I«u)b. ■: 




Ci • ( Ul M 11 / 8lfl2 ): 






tl . (m) lir/aia ); 




1 


raum ( c ); 




x cea<:»> 



E)~6coap~ 



E>-doubJa 

t 



bcdMa, 
a. b; 



lo 



X. XX, Zti 



% • <)oe*)b.x • 0<ac)b.x ♦ (loa<)b,y « ()o&f)b.r: 
X !• 1024; 

If ( t - 0 )' fprlntf {»td«rr»*?B DM<»» by 0 •rrar is bexH* fO: 
zx - <loot}«.x • U<mc)b.x * dene)*.' • (looxJcr: 
it - (ions)cy " (loax)b.i - Uoo«} a, i • (loo«>b.y; 
cx • < lot ) «« / x)««>: 
tr • < Ut ) «t> / *>•«; 

rataro ( c ): 



ml • ) 



>U x.y.x; 



t - (toibl«)c x; 
r • (double)*.*; 

x • wif «•!♦»•»> / «19Z.O: 
ratorn < X ): 



-400- 



fi,1 ¥1-177/161 (7) 



) 

tamp 
lit 

I 



-boc»p to6tt**(». b) 



bi * (tat) ( « • ii 
C| - (1D1) ( fa • 8182.0 ): 
rmim ( c }; 



for< t ' 
I 



<t).i /• I: 



5) -»o id bcaon2{t. a. i 



p • a: 

miIim ■ 0; 

for < I • 0 ; I < • ; ) 

If i «faa((*p).s) > Milu*) hum • *«(«»p>.»): 
If ( «tM({«p}.y) > utlM) »<1M • •b*:(«p).rh 

I 

oobh ■ i: 

•fail* < mi t bus > 4)82) 

t 

•uIm /• X; 
Mp 2: 
cam «- t; 

I 

tf < ommx - I > r*mro: 
P - »: 

for ( I • 0 : t < o ; I** ) 
1 

f«p). i /• oowi: 
(«p).f /• COM4: 

» 



■ •■vS.lp.l J.fe.«. !•.]•!: 



((•L)po»C2.», I); 
!•»/!: 



fp-far (t-O.i-C; fCa-t i 
I 

lf( K»J > 
t 

i • *01; 
tit] ■ *Uh 
iCJ] • t; 

I 

k -wb 
•siU ( k<-i ) 

I 

I - U 
k h *: 

1 

J «- »; 



3- for < ; m<-l : •*+ ) 



I* • (l»t)pa>(lC, (doubles) j 

t.i • i« / t: 

. * toboo^Cl. 0, 0.0)| 

c - cn( I PI / (dowbU) Ul >i 

• * ( <4m*U) f ) - .!.( |_PI / (dubU) Ul ): 

• • tobccop( e, • )t 
for U-«: J<!«/I :J»*) 
I 

for < l»J: t < ■ : !*-!•) 
i 

IP • I ♦ 

tClp] - bcnb(t(l}.D; 
»tiJ • boddftCU.Ot 

I 

• ■ bowl{ m. m ): 

I 



Ibi 



ulpfoid) 

di. i. r; 

I. Mi 
tflttl): 

4i * L0 • * PI / 1024.0; 
i • 0: 

for( i-o : i < 1024 ; )h ) 



*0) • 



■u(i).o.o ) t 



prUtff* TIM oqmIb H">i 
lor( 1*0 : I < ■ j i** ) 



i • ( <dwbl*)i(l].i > / 8162.0: 
oriPtf( * tlva).* • llf i«\ I. i ): 



PT tlM()t 

d hi. i. 10, -I ); 
pr UmO; 

prUur FmpMoep am* In »o*); 

for( 1*0 ; I < ft : ) 

I 

i - ( (doabMt(i).i ) / B1M.0; 
t • < (doabj •)&£!].; ) / lltt.0: 

prlatf( * t[XJH).i • llf. xfttdhf • llf I. l. I, F>> 

I 
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u 








| do«bU raj 












| eoapiai: 






cotplat ewi<H).ci.b().ca»M}.cdt»{) .eoaJO.c.ieO.cpo 


»(). 




cr*() .locoaplaiO ,el«p<) . 






doabta bs ( } *c • r a( ) ~ t 




I> — 


cnapfai Cinti, oj 






co»f>l«i a, bi 

f 






eoaplas »; 






». ra-a. ra ♦ b. raj a. la ♦ b. I»; 






ratara t; 

) 




3>- 


cai>lti caab(a.b) 






coaplai a. b; 






r 

cowlti i; 






». ra»». r* - b. r»: a. la-a. I* - b. f •: 






ratara s; 

1 






eoiplii caal (»,b) 






coaplat a. b: 






f 

coaplai «: 






■. r»*a. ra«b. ra-a. Ia»b. la; 






t. la-a. ra-t>. ta»a. la'b. rat 






ratara a; 

J 





&-coaptai edlv{».b) 
cm>Ui a, b: 

coaplat i: 
doabta a: 



3-doobU c»bi{a) 
coaptat 



4-b. ra » b. U-b. la; 
if { *»*0 ) 

pal at f r\T\T\7«lH«a« bp laro arrorlM\a*): 
a • TUT; 

I 

t. r*«{». r»*b. rt*i, la»b. la)/d; 
i. Ia»(a. ia»b. ra-a. ra«b. Ia)/d; 
ratara i; 



ratara (tqn(«. raaa. ra«a. !■**, la) I : 



doabl* carf(a) 
coapUi a: 

I 

if < (a 
I 



prlaif("\T\7\? art aaro •rror 1 1 1\d') t 
ratara 10.0; 



•la* ratara (ataaZU. la, a. eo))i 



coaplat coaj(a) 
cmjUi a; 



t. ra-». ra: a. Ia*-a. la: 
ratara at 



3)-co«ptei toeoatai(i.r) 
4oabl* »,»: 
1 

utM*> a: 
». ra-aj *.!*•?: 
ratara a; 



i>-i«t cfM i. I. h 

eoapiai 
uuabia f: 

Im I: 

lai a, a*t, la. I. J. ». a. la. lal: 
coaptat t. a, ai 
doabla a. e : 



l>-for (1-0. J-Q: Ka-l: !♦♦) 



/• FFT by aitai eoapiai irpa 
/• adaptad bf S. lolki froa laf(3) 

f* f— 1.0... fft 
/• 1.0.. . rfi 



a - ipoaZ(l); avt-a/t ; 

/« bit ravaraal oparailoa •/ 

If ( K-J ) 

i-aCl]: a(l]*t(J]t «CJJ-t: 



I 

«-a»Ij 

•fella ( t<*J ) 

J— it t/-J; 

) 

j»-a: 



I 

/..•*• Mia loop •••••/ 

tb-for { a-l{ a<»l: •♦♦> 

I 

la 



IpoaZ(a); lal-la/Z: 
a • iocoaplaid.0. 0.0); 
e • coa( •!/(doabla)lal )j 
» • f*ala(ft/(doabU)lal); 
* - toeoapl*i( e. a ) ; 
for U-0i JCta/t; !♦♦> 

for ri-Js Kai I*-!*) 

lp-1+lali 

l .ewlta[l*].o>: tflal-eaaMttll.O; 
Btl)*c-sa(stll.i)t 

1 

a*caal( a. a); 



TB e ten 



llacl&da <atdto.b> 
tlacloda <aatb. k> 
lincloda <coapla>. b> 

©-void effi(f.l,f) 
cuopl as "tj 
i.t l.f; 
I 

lat a, a*; 
coaplai l, a. 



a * (ta«)ao*(2.0.(doubla)l): 
a*t • a / I; 

Ip-for (1-0. J-0: Ka-l : !♦♦) 
I 

IM U-J ) 

1 t • ifilt 

afi) • x(J): 
ttJi • w 

I 

a • a*ti 
ablla ( k<-J ) 
I 

i - t: 
k h U 

I 

j ♦- b: 

1 



w t ted <»r 
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fo -tor 



T" 



l.l •!•/*$ 

■ • taco«»U»(i.O, 9.(1): 

c - coj< an/ (tfoi>bi*> i«t ); . 

, • ( (tf M Sl«) f ) • tJo( a.Pl / (tosfcU) Ul )l 
« • toeo*9l*i( c, • )i 
for <J*0; #<!•/■ : J**) 

for ( l-J: I < i : t U ) 
I 

IP - I ♦ Ul; 

x(ip) • Ct»b(l£i].Oi 

s[U - eiM(i(t],tl; 

» 

• • cibK >. ■ ): 

I 



( f < 0 > for { J • 0 : I < » : ) 



t(i). ■ /- (do^bU) •: 



o-Z* .no Z-a/Z 




j-j+i.i-ui 



% 8 ® 
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